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ABSTRACT 


A circular cylinder immersed in an infinite fluid and oscillated 
perpendicular to its axis is acted upon by a fluid dynamic force. This 
force is generally considered to be composed of two components, a drag 
and an added mass component. The relative importance of these are 
functions of the frequency and displacement of the cylinder. 

It was the purpose of this study to determine experimentally these 


two components as a function of the ratio of the displacement amplitude 


to the cylinder diameter, and the Reynolds number based on the maximum 


velocity. 
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I. INTRODUCTION 


With the increased deployment of submerged ocean structures the 
need to predict the forces exerted on bodies immersed in an oscillating 
flow is of practical importance. Studies conducted by Reid and 
Bretschneider [1] on submerged structures stress the necessity of 
having these forces forecast, given a stated set of wave parameters. 

A number of studies have been carried out to determine the values 
of the added mass and drag coefficient for a pile in ocean waves. Many 
of these are represented in a correlation of drag coefficient versus 
Reynolds number presented by Wiegel [2]. However, his plot shows an 
order of magnitude scatter in the values. 

EM basic studies of forces exerted on bodies in oscillating flow 
have been conducted. Of these, however, the most well known was 
carried out by Keulegan and Carpenter [3]. They conducted tests with 
a stationary cylinder placed at the node of a standing wave and obtained 
a correlation of Cm (mass coefficient) and Cd (drag coefficient) versus 
Um T/D, where Um denotes the maximum velocity, T the period, and D the 
cylinder diameter. It was suggested by Keulegan and Carpenter that the 
parameter, Um T/D, or equivalently, the relative fluid displacement, A/D, 
was the parameter of primary importance; the effect of Reynolds number 
was assumed to be of little ae on the values of the drag and added 
mass coefficients. It is noted, ee that the test method used by 
Keulegan and Carpenter allowed no control over the Reynolds number so that 


they had no convenient method of testing this influence. 





Heinzer [4] conducted visual studies of an oscillating circular 
cylinder and concluded that the wake characteristics were primarily 
dependent upon amplitude of oscillation vice Reynolds number. 

Sarpkaya and Garrison [5] investigated unidirectional flow past a 
cylinder with constant acceleration. Their work lead to the fact that 
drag and inertial forces are representable as functions of the relative 
displacement of the fluid. These studies were conducted in the range 
of Reynolds numbers near 10° and no attempt was made to investigate 
the effect of the Reynolds number. 

In the present study various circular cylinders were oscillated with 
simple harmonic motion in water at rest. Both the amplitude and the 
frequency of oscillation were controllable so that both the effect of 
Reynolds number and relative displacement of the fluid could be studied. 
However, time limitations restricted the experimental program to a 
series of runs carried out at one fixed Reynolds number with increasing 
relative displacement. Also, at a fixed value of relative displacement 


the effect of Reynolds number was determined. 





ТТ. METHOD OF APPROACH 


This: study utilized test cylinders of 4 inches, 2.75 inches and 
1.5 inches diameter, 16.5 inches long mounted on a dynamometer. The 
dynamometer was baste composed of a cantilever beam affixed with 
strain gauges. The strain gauge bridge error was read out on a strip 
chart recorder calibrated to read directly in pounds force. 

The force traces were read visually at twenty intervals over a 
complete cycle and this numerical force data used as an input ES a data 
reduction program (Appendix A and B) based on Morison's equation. 

Morison et. al. [6] consolidated the three part equation of McNown 


and Wolf [7]. 


Бе Ао р ЖЫ) +$ Px ds +Z Cd D p U|v| (1) 


by combining the first two terms to yield 


du di 
Е = Сшр 49 #5 29 Ом ПШ (2) 


where F --- Force per Unit Length 

Ao --- Circular Cross Section Area of Cylinder 
D --- Body Dimension Normal to Flow 

р === Density of Fluid 

Е --- Time 

Са --- Coefficient of Drag 


Cm --- Added Mass (Inertia) Coefficient 


U --- Velocity at Time t 





Equation (2) is assumed valid and is used as the basis of the data 
reduction. 
In the present instance the velocity of the cylinder is represented 


by: 
U = -Um cos Oo t (3) 


where Um denotes the maximum velocity, T the period, and а + 21/Т. 
The total force acting on the cylinder per unit length is in general 


, 


given by 

EE AD O) (4) 
Grouping the variables on the basis of dimentional reasoning and 
introducing a = 21t/T gives 


F UmT UmD 
z =f, De ee (5) 
Um D 





UmT 


where E is a Reynolds number and ——, which can also be 





expressed as will be termed the relative displacement. 


277A 
D э 


Because of flow symmetry and the periodic nature of the force 
СЕБЕ вл) (6) 


it is possible to express the force coefficient in a Fourier series, 


- 


—— = A. sin a + A. sin 3a + A. Sin 5a +... 


1 3 5 


(7) 


+ В, cos a + В, cos 3a + В. cos 5a +... 


10 






where the coefficients Ay» By are independent of a and at most 


functions of Reynolds number and relative displacement. Fourier analysis 


may be used to determine the coefficients as: 


| 27 Е 
a 225 тэш da (8) 
0 oUm D 
and 
21 
в, == F => Na da (9) 
0 pUm D 


Once obtained, the dependence of these coefficients on Reynolds number 
and relative displacement may be established provided the data are 
sufficient. 

The general formulation, equation (7), may be reconciled with 
Morison's equation (2). Introducing U from equation (3) into 


equation (2) 





5 =F om + 2% sin a - 55 cos a cos a (10) 
pUm D i 


By the rule of Fourier 


P 
6 | cos a | cos a cos Nada 


27 2 
| cos Nada 
6 


Ms 


COS O | COS Q 


Z, 
li 
© 


7 89 + а, cos a + а, cos 2a + a, cos 3a + ... 


where 


(11) 


11 





The first three non-zero coefficients would then be: 


a, = 8/37, a4 7 8/157, а; = 8/105т (12) 


Introducing equation (12) into equation (7) with 


1 1/31 


е3 
| 


De Do аз/ а, [Bj] 
; (13) 





Be = В, ~ a,/a, |81) 
yields 
F я 5 4 
= A. sin a ҒА, sin За + A. Sin 50 +... 
2 1 3 5 
pUm D 
(14) 
t 
E t t 
+ Bi | cos a | cos a + В) сов Dot B¿ cos 5@ +... 
Equations (14) and (7) may be compared. Writing 
T MEDIE sin 3a sin 5a 
4 u Um = x Az sin a 5 sin a d (15) 
and 
? ` 1 
Cd : B4 cos 3a Be COS a 
— = -B == = + °°. (16) 
2 1 | cos a | cos a | cos a | cos a 


Thus, if the coefficients Аз» Ac and В, Б В. vanish, the 


values of Cm and Cd remain constant for all phases of the cylinder 


motion 


_ 2 Unt 
Cm = E BUS А, (17) 





Substituting A, from equation (8) yields 


1 


UmT E F sin a da 


2 

Cm = — (18) 

т^ : 0 pUm^D 

, 
са -72B, (19) 
substituting from equations (9) and (13) yields 
| 21 : 
Cd - = Е sin a da (20) 
0 pUm D 


If the coefficients do vary with the phase a, the values given by 
equations (18) and (20) are weighted averages. With this possibility 


in mind it is preferable to adopt 


F , 








z = 4sin a * Bj cos a | cos a | + error (21) 
pUm D 
or 
= = 2 Сш E а - 54 | cos a | cos a + error (22) 
pUm D Um 
, 
where А,, B4 » Cm and Cd are constant and error has the value 


error = A, sin 3a + A, sin 5a 


~ 


+ B, cos 3a + В. cos 5a 


This error is obtained by subtracting the computed value of Ay sin a 


! 
апа В, | cos a | cos a from the observed F/pUn D. 


13 





The foregoing formulation was developed by Keulegan and Carpenter 
[3] and equations (18), (20) and (22) were the basis for the computer 
data reduction program developed in Appendix A and presented in 


Appendix B. 
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ТТТ. APPARATUS 


A. TEST BASIN 

A test basin, 16 feet long, 4 inches deep, and 16.5 inches wide was 
constructed from 3/4 inch AC marine plywood. The side walls were 
supported at approximately 20 inch intervals with 1.5 inch thick plywood 
gussets and two 2 by 6 inch pine stringers along the top edge of the 
tank, Figure l. 

The presence of considerable width-tolerance variance when the tank 
was filled, necessitated the addition of three 2 by 1 inch channel 
braces, welded in an inverted U-shape, placed over the top of the channel. 
The brace design allowed for unobstructed motion of the dynamometer 
along the length of the tank, an advantage not afforded by structural 
cross braces along the tank top, Figure 2. 

All joints were sealed with RIV. and inside surfaces were treated 
with five coats of penetrate and five coats of epoxy paint. Two 3/4 
inch thick plexiglass windows, Figure 3, (20 inches by 12 inches) were 


installed near the tank midsection. 


B. DRIVING MECHANISM 

The driving mechanism consisted of a "Dynamatic" electric motor with 
an eddy-current clutch and Mark III controller. This system utilized a 
closed loop feedback to linearize and stabilize its shaft speed. The 


basic control system used is presented schematically in Figure 4. 
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Figure 4: Feedback Loop 


More specifically, the controlled quantity is the shaft speed, the 
reference input is a potentiometer and the error detector is a summing 
amplifier. The controller is a combination of a power amplifier, 
motor and eddy-current clutch. The measure of shaft velocity is 
accomplished by means of a tachometer generator, producing an output 
voltage proportional to shaft velocity. 

The summing amplifier detects differences between the desired speed 
and actual speed. This med Wa fed to the power amplifier which 
modulates the clutch coil current. If the speed is too low, the coil 
current is increased, causing the eddy-current clutch to transmit more 
torque to the load. The result is an acceleration of the shaft to 
red speed. 

The above speed regulation operates for positive loading only. 
Since the minimum coil current was zero, negative loads caused free 
rotation of the output shaft. 

Both positive and negative loads were experienced by the output 
shaft due to the connecting rods weight and carriage inertia. This 


lead to uneven output shaft speed, particularly at low RPM, and 


18 





necessitated the application of a continuous positive load. This load 
was supplied through a 3 inch wide 17 inch diameter brake drum, affixed 
to the face plate, with a 120° leather lined brake shoe riding on it. 
The assembled motor, face plate and brake were mounted on a separate 
platform, Figure 5, to isolate any inherent mechanical vibrations from 
the tank structure. A 7 foot 8 inch long by 1.5 inch diameter thin 
walled aluminum tube was then used to connect the driving mechanism to 
the carriage. This long driving rod produced a near sinusoidal motion 


of the carriage. 


C. CARRIAGE 

The carriage as shown in Figure 6 provided a movable mounting plat- 
form which was positioned over the test tank. 

The carriage consisted of a U-shaped, welded aluminum structure 
fitted with four 3/4 inch linear bearings located at the ends of the 
hollow cylindrical side members. The sides of the platform were 
constructed of two inch aluminum pipe fitted with aluminum bearing blocks 
at each end. The structural transverse member of the carriage was a 2 
by 3 inch aluminum box section welded to the end of each pipe. At the 
opposite ends of the pipes the bearing blocks were connected by a 3/8 
inch diameter transverse shaft to accommodate the vertical struts which 


~ 


Supported the cylinder. I 
The U-shaped platform was mounted through linear bearings on two 
case hardened 3/4 inch shafts. The shaft ends were supported by two 


2 by 6 inch steel channels that transversed the width of the tank. These 


channels were clamped to the stringers along the top of the side walls. 
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The movable part of the carriage consisted of two long arms pivoted 
near their centers and mounted on a drum as shown in Figure 7. These 
arms were designed to hold the test cylinders at a maximum centered depth 
of 18 inches submergence. The leading edges of the tapered 1/16 inch 
thick polished support arms were sharpened to reduce drag and experimental 
errors. 

A counter balance was placed on the support arms opposite the test 
cylinder and adjusted so that the center of gravity of the system was at 
the pivot point. As a result of this balance adjustment the only force 
(or moment) acting on the arms was due to the fluid; all other inertial 
loads were balanced out. 

A third arm, half the length and perpendicular to the supporting 
arms, was welded to the center drum. This arm served as a 2:1 mechanical 
amplifier of force and converted the torque about the rotational center 
of the drum to a linear force. This force acted, through a pinned ball 
bearing linkage to a small cantilever beam equipped with strain gages, 


Figure 8. 


D. CANTILEVER BEAMS 
Two cantilever beams were constructed, Figure 9; 


1. 0.25 inch wide by 0.35 inch high with action arms of 
2.5 inch and 4.5 inch 


2. 0.6 inch wide by 0.45 inch high with action arms 
В ımehwors 3.5 inch 


wo attachment points were provided on each beam so that a large range 
of forces could be measured without overstressing the beam. Stiffness, 
and/or sensitivity, could be adjusted to reduce natural frequency vibra- 


tory motion by shortening the attachment point or by a change of beams. 
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Four strain gauges were applied to each beam to provide a means of 
measuring forces imposed upon the cylinder and transmitted mechanically 


through the balance assembly to the beam. 


E. TEST CYLINDERS 

Three cylinders as shown in Figure 10 were constructed to provide 
maximum flexibility in testing over a wide range of frequency and 
displacement. All cylinders were fabricated from plexiglass and slotted 
to accommodate one end of the vertical support arms. Bulkheads were 
placed adjacent to the slots to provide attachment points for the ends 
of the arms. Silicone sealent was used to prevent leakage after the 
arms were attached. This allowed the cylinder to be unflooded resulting 
in a minimum fixed mass. End pieces as shown in Figure 11 were also 
constructed of plexiglass. The ends were fitted into the cylinder with 
1/8 inch "0" rings to provide the necessary seal. The outer surfaces 
of the end pieces had a labyrinth seal cut into the face and 1/16 inch 
plexiglass end plates, one inch larger than the cylinder. The labyrinth 
seal and end plates were provided for purposes of reducing the leakage 


around the ends of the cylinder. 
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IV. TEST PROCEDURE 


A. SYSTEM SELECTION 

Primary consideration was given to the selection of the cantilever 
beam to be used with each of the test cylinders, trade offs being stiff- 
mess and sensitivity. It was concluded that the 4 inch and 2.75 inch 
diameter cylinders would be run on both beams and the 1.5 inch diameter 
cylinder tested only on the more sensitive of the two beams. This 
procedure yielded force traces that could be cross correlated for the 
4 inch and 2.75 inch cylinders with the better trace utilized for data. 

Once components were selected, they were assembled, balanced, aligned 
to the tank axis and clamped to the tank edge. The system was then finely 
balanced in the dynamic mode by adjusting the lever arm of the counter 
balance weight after observing strip chart traces of high frequency 
oscillations in air. The counter balance weight was adjusted so that no 
force was observed when oscillating the carriage in air. The tank was 
then filled to a level 4 to 6 diameters above the test cylinder and 
the desired amplitude of oscillation and frequency set. Amplitude 
settings, i.e., connecting rod eccentricity on the face plate, were in 
half diameter increments. The frequency setting was either random when 
observing the Reynolds number parameter or adjusted to maintain a given 


Reynolds number for studies of the amplitude/diameter parameter. 


B. SYSTEM CALIBRATION 
1. Force Calibration 
A length of fine nicrome wire was attached to the test cylinder, 


fare lead over a pulley arrangement, attached to the top of the tank. 
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After a half hour warm up period, the four arm strain gage bridge was 
electrically balanced in an unloaded state. Then accurate weights were 
hung from the wire and the resulting deflection of the recorder "T noted. 
The strain gauge bridge and recorder were calibrated in both the positive 
and negative direction. The resultant deflections of the recorder pen 
were accurately plotted against the force imposed, Figure 12. The force 
calibration proved to be very linear. 

This calibration was conducted prior to all test runs. 

2% TER Marker Calibration 

Timing of the moment of maximum velocity was determined by 

timing contacts on the drive wheel. The ability of the event marker to 
record the frequency of the imposed motion and the phase angle between 
the force and motion was dependent upon the accuracy of the chart speed. 
However, this speed was checked against internal and external timers and 


no discrepancy was found. 


C. DATA RUN 

The system as selected was set in motion and recorder started. To 
avoid wave reflections from the ends of the tank, the system was 
oscillated for less than ten cycles and secured. No appreciable standing 
wave was generated under this procedure and the system was allowed to 
remain at rest for fifteen minutes before a new data set was initiated. 


After each data run, the force calibration and event marker position 


were checked. 
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D. STRIP CHART READ OUT 

A mean curve was fared through the force recording to smooth out 
natural frequency fluctuations and bearing noise. This curve was then 
divided into twenty divisions, Figure 13, and mid-division force values 
were read and correlated to its respective t/T for inputs to a data 


reduction computer program, Appendix B. 





Figure 13: Sample Force Read Out 
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V. RESULTS 


The purpose of the experimental program was to study the effect of 
the parameter A/D on the added mass and drag coefficients at a fixed 
Reynolds number. Also a limited number of runs were conducted to deter- 
mine the effect of Reynolds number on these coefficients at a fixed value 
of relative displacement, A/D. This latter series of runs was particularly 
significant since in the past it has been taken for granted that little 
or no Reynolds number effect occurs in flows of the present type. 

КЫ йге 14b and 15b show the total calculated force, based оп 
Figures 14a and 15a, and the observed force for two typical test runs. 
These figures show that the representation of the total force as the sum 
of the drag and inertia components is adequate for practical application. 
Figure 14 corresponds to a case dominated by inertia while Figure 15 
corresponds to a larger value of A/D where the drag component is the 
greatest. The error computed in Appendix A is the difference in the 
observed and calculated values plotted in Figures 14 and 15. 

Figures 16 and 17 show the results for Cd and Cm plotted as a 
function of A/D for a smooth and een cylinder. Also, Keulegan and 
Carpenter's mean value line [3] is plotted for comparison. Keulegan 
and Carpenter conducted their studies with stationary cylinders and 
standing waves and, therefore, had no independent means of controlling 
Reynolds number for a fixed A/D. In general, however, their studies 
corresponded to fairly low Reynolds number. 

It is noted that the mean Cm line taken from their work has been 
reduced in value by 1.0 to account for differences in test procedure 


(stationary fluid versus moving fluid). 
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Figure 14b: Calculated vs Observed Force 
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It can be seen that Cd values which correspond to rather high Reynolds 
number, fall below Keulegan and Carpenter's results. However, the series 
runs which were made at fixed A/D with increasing period (decreasing 
Re ds number) show that the values tend toward the low Reynolds number 
results of Keulegan and Carpenter. It is surprising, however, that the 
effect of roughness is to increase the drag coefficient so as to approach 
Keulegan and Carpenter's mean line. 

Figure 17 shows the variation of Cm with A/D along with the results 
of Keulegan and Carpenter. It is noted that the results for rough 
cylinders and for higher Reynolds numbers encountered in this study fall 
above those of Keulegan and Carpenter. It is interesting that roughness 
tends to have little effect on Cm in the present study. 

The results for Cm and Cd plotted against A/D is presented in 
Figures 18 and 19, respectively. No effort was made in these plots to 
keep the Reynolds number constant. As a result the Cd plot shows a 
large amount of scatter; the Cm plot, however, appears to correlate 
fairly well without regard for the Reynolds number. 

Figure 20 and 21 show the added mass and drag coefficients plotted 
against A/D for a fixed Reynolds number. It is noted that the scatter 
caused by the non-constant Reynolds number has been suppressed, and the 
results correlate well. 

Keulegan and Carpenter's mean lines also have been plotted for 
comparison on Figures 20 and 21. The magnitudes of individual points 
do differ but the same trends are prevalent throughout. This difference 
is attributed to experimental differences (particularly Reynolds number) 
between studies and inherent experimental error. The critical region of 


A/D remains in the region of 2.25 to 2.5. 
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Figures 22 and 23 show the effect of Reynolds number on the added 
mass and drag coefficients, respectively. Figure 21 shows that the added 
mass coefficient is almost independent of Reynolds number and highly 
dependent upon A/D in the range tested. However, Figure 23 shows that 
Cd is strongly affected by the Reynolds number and generally it shows a 
decrease with increasing Reynolds number. 

Keulegan and Carpenter [3] state that apparently no correlation 
exists between the two coefficients and Reynolds number. However, on 
plotting their data in the format of Figures 22 and 23, the same general 
trends are observed; Cm appears to be independent of Reynolds number 
while Cd is strongly dependent. This data is plotted in Figures 24 and 
25. It appears that Keulegan and Carpenter overlooked this rather 


Significant feature. 
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VI. CONCLUSIONS 
Total force acting on a cylinder oscillating harmonically 
is comprised of two components, drag and inertia. 


The total force may be mathematically modeled using 
coefficients of added mass and drag. 


Cm is primarily a function of relative displacement A/D 
and appears to be nearly independent of Reynolds number. 


Cd is a function of both relative displacement and Reynolds 
number. 
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APPENDIX A 


A. FORMULATION OF THE COMPUTER PROGRAM FOR DATA REDUCTION 

A computer program was used to reduce the data. The purpose of the 
program was to ТЕКС Cm, the inertia coefficient, and Cd, the drag 
coefficient. The starting point for the problem was Keulegan and 


Carpenter's [3] formulation for these parameters. Specifically: 


I 57 Е зіп F sin a 
Cm = 


d a (1A) 
Е D O Aline D 
and 
Cd = -3/4 ç EA da (2A) 
0 ойт D 


The above integrals were transformed into. the following forms: 


2T 2Tt 

Cn = —<—— F sin (—) d(t/T) (3A) 
DT p A 3 = 

pt 
-3T 2T £ 

Са -------- F cos С--) 4 («/Тт) (4A) 
8: p Duc 2 £ 


where A denotes the amplitude of the oscillatory motion and d(t/T) 
was conveniently chosen as 0.05. 
The corresponding error for each data set was computed by the 


following relationship: 
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error 





[> E ТН э d S (225) 


Т 6Tt 
+ Е сов (----) J(t/T) 
2510 Dacus 2. Т 


(5А) 


0.02 2 Ж (225) (ey) ав (Сот, 


where: 
T = period in seconds 
D = cylinder diameter in feet 
L = cylinder diameter in feet 
t = time 


This error function is obtained by subtracting the computed values from 


the observed. 
The program employed time increments of 1/20 T. where t/T takes on 
values from zero to one. For each run at constant period it was possible 


to fix t/T increments for each data set. 
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B. FLOW CHART OF COMPUTER PROGRAM 


DEFINE 
PARAMETERS 


INITIALIZE 
CONSTANTS 


INITIATE 


DATA 
SETS 





READ 
PARAMETERS 






COMPUTE 
“| COEFFICIENTS 
i БАЗЕ 


1725 


COMMENCE 
DATA 
REDUCTION 


COMPUTE 
TRIG 
ARGUMENT 
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COMPUTE 
FORCE 


PRENT 
COMPONENT 


COMPUTE 
Cd Cm 


COMPUTE 
ERROR 





STOP 
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APPENDIX C 


A. CALCULATION OF INERTIAL AND DRAG FORCES 


Given 
T, Cd, Cm, D, Am, and 
p = 1.93 slugs 
F = F. + F = Cn o Ao 0 + 2 сар у |u| 
I D Р dt 2 0 
K аа u 
I dt á dt 
U = -Um cos о Е ші. Um o sinc t 
Ш - 21/1 
Ат 
Um = 2 7 = 
Therefore: 
ош, zm Am SADO owt 
dt T Т 
An? И 
= 5 sin co t 
T 
Therefore: 
F_ = Cn оп p? 4n? Am en 2Tt 
I 4 T ab 
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Let ae 
e а T a(t) 
with | О-о ЕТ 
2 
Спо т р Am . 
Е ------------ спа 
I pe 


F = баре и [ul 


substituting for U 


E = = Cd Do Um” | cos o t | cos o t 


D 


substituting for Um 


e 


EB = = Ms D p eh cos ot | cos o t | 
D 2 ge 
therefore 
2 m 
Fy = - Cd = = L cos а | cos а | 


when a is substituted for ot. 
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T = 1.60 seconds 


p = 1.93 slugs 


= 0.508 sin a 
= -0.110 cos a 
— 5 
Time Imertial 
Period Force 
0.25 0.979 
0.75 0.230 
1.25 0.360 
1.75 0.453 
2425 0.502 
2.75 0.507 
3.25 0.453 
3.75 0.360 
4.25 0.230 
ATS 0.080 
5425 -0.080 
5.75 -0.230 
6.25 -0.360 
6.75 -0.453 
7225 -0.502 
7. 75 -0.502 
8.25 -0.453 
8.75 -0.360 
9.25 -0.230 
9.75 -0.079 


CALCULATED DRAG AND INERTIA FORCE FOR: 


— m => ——— — . — — = — F  - — = | емел эч 


COS а 
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Cm = 1.1722 
A = 0.1667 ft. 
Drag Calculated 
Force Force 
-0.107 -0.027 
-0.087 0.143 
-0.055 0.304 
-02029 0.430 
-0.003 0.499 
0.003 0.505 
0.028 0.475 
0.550 0.414 
0.087 0.318 
0.107 0.187 
0.107 0.048 
0.087 -0.143 
0.055 -0.304 
0.023 -0.430 
0.003 -0.499 
-0.003 -0.505 
-0.023 -0.475 
-0.055 -0.414 
-0.087 -0.318 
-0.107 -0.187 


Cd 


0.7947 


0.3333 ft. 


Observed 


Force 


0.05 
0.20 
0.33 
0.43 
0.50 
0.50 
0.50 
0.45 
0.35 
0.19 
0.05 

-0.10 

-0.25 

-0.40 

-0.45 

-0.50 

-0.50 

-0.45 

-0.35 

-0.10 





CALCULATED DRAG AND INERTIA FORCES FOR: 


——— с | SS ee 


T = 1.39 seconds Cm = 0.9838 Cd = 1.8856 
p = 1.93 slugs Am = 0.3437 ft. D “= 0,299) ie, 
Fr - 0.667 sin a 
Fo = -3.17 cos @ | cos au 
Fo = Fr ЕЕ Fy 
Time Inertial Drag Calculated Observed 
Берсе Force Force. Force Force 
0.025 0.087 -0.983 -0.896 -0.900 
0.075 0.252 -0.800 -0.548 -0.600 
0.125 0.393 -0.504 -0.111 -0.300 
0.175 0.496 =0. 207 07289 0.000 
0.225 0.550 -0.024 0.525 0.300 
02275 0.550 0.024 0.574 0.600 
0.325 0.496 0.207 0.703 0.900 
0:375 0.393 0.504 0.897 1.200 
0.425 09292 0.800 1.052 0.900 
0.475 0.087 05983 1.070 0.700 
207525 -0.087 0.983 0.896 0.900 
0:575 -0.252 0.800 0.548 0.700 
0.625 00393 0.504 0.111 0.250 
0.675 -0.496 0.207 -0.289 -0.250 
0.725 ae 0.024 -0 525 -0.600 
02775 -0.550 -0.024 -0,574 -0.800 
0.825 -0.496 -0.207 -0.703 -1.100 
0.875 -0.393 -0.504 -0.897 -1.100 
0.925 -0.252 -0.800 zm -0.700 
09.75 -0.087 -0.983 -1.070 -0.650 
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DATA RUNS FOR REYNOLDS NUMBER 2.7 x 10 


APPENDIX D 


A ems rm ee eee 


R4 


R6 
К/ 
R8 
R9 
R10 


R12 
R13 
R14 
R15 
R16 


R18 
R19 


R21 


Cyl. Dia. 


о. 
1.5 
1.5 
1.5 
1.5 
15 
1.5 
1.5 
2275 
2275 
2275 
275 
21515 
2.75 
DU 
2.75 
2:15 
4.0 
4.0 
4.0 
4.0 


A/D 


I 
2.0 
2,5 
3.0 
35 
4.0 
4.5 
5.0 


N N 9 Fa PF Fa ы-і О 
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4 


0.46 
0.60 
0.78 
0.91 
1.05 
1221 
1.37 
1.56 
0.75 
0.75 
1.0 

0.96 
1.50 
1.45 
2.06 
253 
2.57 
0.53 
1.50 
1.06 
2.15 


Cm 


1.4692 
0.2134 


-0.1047 


0.3988 
0.2462 
027225 
0.0249 


x. 35 39 


1.4791 
1.4594 
1.2939 
1.2832 
0.8159 
0.8762 


-0.6060 


0.1127 


-0.0085 


1.4959 
1.5350 
1.4632 
1.3128 


Cd 


1.7197 
2.1625 
2.4826 
2.0726 
2.1140 
2.0300 
2.0825 
2.2148 
PESTO 
1.0918 
1.1664 
0.9221 
1.7550 
1.9681 
2.9893 
23252 
2.4444 
1.6421 
0.9485 
1.1187 
1.1295 





A A O A E A A < = m= Í = E 


Run 


о о Y Оо іл г QQ > ~ 


10. 
11. 
ша: 
E. 
14. 
ШӘ? 
16. 
17. 
18. 
19. 
20. 
21. 
22, 
23. 
24. 
29. 
26. 


A/D T 
172 0.84 
1/2 0.60 
1/2 0.48 
1/2 0.40 
12 0.33 
1 0.83 
"d 0.61 
il 0.48 
1 0.39 
1 0.33 
155 1259 
ives 0.83 
1.5 0.62 
1.5 0.48 
2 1.30 
2 0.85 
2 9 1652 
2.5 0.85 
25 0.55 
1/2 1.26 
1/2 UNOS 
Il 12936 
1 0.86 
1.5 1.35 
ds 1.05 
T 0.84 
Runs 1 - 19 
Runs 20 - 26 


Cm 


1.4201 
1.4358 
1.4395 
1.4244 
1.2812 
122561 
123027 
13375 
1.3154 
1.2295 
0.9838 
0.9190 
1.0179 
0.9710 
0.8083 
0.7454 
0.5224 
0.7891 
0.9003 
1.4857 
1.2976 
1.4547 
1.3894 
1.1263 
1.0685 
1.1361 
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Cd 


1.0556 
1.102 
0.7405 
0.3118 
1.2897 
0.9318 
0.7524 
0.2872 
0.7546 
0.9274 
1.8856 
1.1842 
1.0826 
1.0334 
1.3861 
1.2833 
1916231 
1.1180 
0.8348 
1.2957 
1.6856 
0.8183 
0.8546 
Е. 
1.0063 
0.8950 


2.75 inch diameter cylinder 


4.00 inch diameter cylinder 


Re x 10. 


T 61 
225 
2,02 
3.38 
4.10 
3226 
4.43 
5.63 
6.39 
8.19 
272 
4.89 
6.54 
8.45 
4.16 
6.36 
5 ПЕР 
7.95 
В 
2,26 
3.43 
4.20 
6.64 
6.35 
8.17 
10.10 
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